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The  Effects  of  Space  Charge  in  a  Hypersonic 
Magnetohydrodynamic  Power  Generator 


Rene  J.  Thibodeaux 

Air  Force  Research  Laboratory,  Wright-Patterson  AFB  Dayton,  Ohio,  45433-7251 


This  paper  will  explore  a  new  MHD  generator  design  that  uses  space  charge  to  create  an 
axial-symmetric  electric  field  in  the  region  between  concentrically  cylindrical  channels.  The 
space  charge  fraction  is  shown  to  have  a  significant  effect  in  the  current  density,  power 
density,  electric  field,  output  voltage,  efficiency,  pressure  drop,  enthalpy,  and 
electromagnetic  Reynolds  number  of  the  generator.  The  radial  symmetry  produces  simple 
analytic  solutions  for  these  equations  over  a  wide  range  of  gas  flow  speeds.  Adjustments  in 
the  space  charge  can  be  used  to  regulate  the  output  voltage  of  the  generator  for  a  constant 
magnetic  field.  This  technique  can  eliminate  the  need  for  a  massive  full  output  power 
conditioner.  Space  charging  current  (electron  beams  or  corona  ionization)  is  calculated. 
Also  the  space  charging  energy  is  used  to  plot  the  system  weight  benefit  of  using  space 
charge  voltage  regulation  versus  voltage  regulation  using  full  output  power  conditioning. 
Plots  of  the  generator  output  voltage,  electric  field,  efficiency,  and  electromagnetic  Reynolds 
number  are  also  presented. 


Nomenclature 

J  =  current  density 

U  =  gas  flow  velocity 

E  =  electric  field 

B  =  magnetic  field 

/  =  current 

ct  =  conductivity 

p  =  charge  density 

Ke  =  electron  mobility 

Ki  =  electron  mobility 

K  =  load  factor 

v  =  collision  frequency 

a  =  electron  density  fraction  -  {electrons/cm3}/ {molecules/cm3} 

P  =  positive  ion  density  fraction  -  {positive  ions/cm3}/  {molecules/cm3} 

X  =  negative  ion  density  fraction  -  {negative  ions/cm3}/  {molecules/cm3} 

5  =  space  charge  fraction  (with  positive  ions  only)  5  =  (a  -  P)/a 

A  =  space  charge  fraction  (with  negative  and  positive  ions)  A  =  (a  +  %  -  P)/a 

I.  Introduction 


MAGNETOHYDRODYNAMICS  has  been  well  demonstrated  in  large  scale  supersonic  laboratory 
demonstrations  since  the  1950’s.  Small  scale  shock  tube  and  wind  tunnel  experiments  to  investigate 
hypersonic  MHD  phenomena  has  also  been  the  performed  since  1950’s,  with  the  most  recent  well  known 
work  focused  on  the  Russian  of  the  “Ajax”  vehicle  concept  first  revealed  in  the  1990’s.  The  Air  Force  Research 
Laboratory  (AFRL)/Propulsion  Directorate  initiated  the  Hypersonic  Vehicle  Electric  Power  Systems  (HVEPS) 
program  in  2001  to  investigate  the  realistic  operation  of  MHD  channels  at  hypersonic  gas  flow  speeds  and 
temperatures  equivalent  to  scramjet  engine  conditions  lasting  for  durations  of  several  seconds.  Development  of  this 
technology  could  make  air-breathing  MHD  auxiliary  power  systems  available  for  hypersonic  vehicles  for  a  wide 
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variety  of  high  power  applications.  Results  of  these  tests  would  conclusively  answer  whether  or  not  scramjets  could 
directly  power  MHD  generators.  The  test  rig’s  mass  flow  rate  was  the  known  limiting  factor  in  the  output  power  of 
the  generator.  The  major  unknown  was  the  non-uniformities  in  the  plasma  caused  by  the  hypersonic  gas  flow. 

HVEPS  was  a  collaborative  effort  between  General  Atomics  (GA),  LyTec  LLC,  Pratt  &  Whitney 
Rocketdyne  (PWR),  United  Technologies  Research  Center  (UTRC),  and  NASA.  The  HVEPS  program  studied  the 
feasibility  of  using  MHD  generators  as  electrical  power  sources;  thereby  overcoming  the  limitation  in  scramjets 
which  lack  rotating  shafts  in  the  direct  gas  flow  and  therefore  can  not  directly  power  conventional  rotating 
generators.  The  UTRC  test  was  one  of  the  final  tasks  of  HVEPS  which  has  include  vehicle  system  studies,  magnet 
and  generator  designs.  This  testing  effort  was  performed  at  the  UTRC  Jet  Burner  Test  Facility,  Hartford, 
Connecticut.  The  test  setup  consisted  of  the  scramjet  combustor  coupled  with  the  MHD  generator  in  the  exhaust 
stream.  The  MHD  generator  was  a  constant  cross-section  rectangular  geometry  with  a  conventional  diagonal 
electrode  lined  channel.  For  this  demonstration,  the  system  used  a  helium  cooled  superconducting  magnet  (3  Tesla 
capability)  supplied  by  NASA/Marshall  Space  Flight  Center.  The  combustion  fuel  was  ethylene  mixed  with  NaK 
seeding  and  oxygen  to  raise  the  combustion  temperature  above  4000  Rankine  for  sufficient  electrical  conductivity. 

The  goal  of  the  tests  was  simply  to  measure  the  voltage/current  during  a  4  second  burst  created  when  the 
seeding  released.  Power  was  produced  immediately  every  time  seeding  was  applied  to  the  combustor  in  6  test  runs 
made  over  the  two 
separate  days  of  testing, 
giving  a  high 

confidence  that  the  test 
results  are  reproducible. 

The  test  results  were 
used  to  map  the  power 
curves  for  various 
electrical  load  levels. 

The  voltages  varied 
from  100-400  VDC  and 
the  power  levels  ranged 
from  2-15  kW 

depending  on  loading. 

The  real  power  levels 
were  limited  by  the  test 
rig  mass  flow  rate, 
instantaneous  electrical 
conductivity  of  gas,  and 
the  magnetic  field. 

Figure  1  HVEPS  Scramjet  Driven  MHD  Generator  with  Helium  Cooled  Magnet 


Figure  2  Voltage  and  Current  Measurement  from  HVEPS  Scramjet/MHD  Generator 


The  system  benefit  of  applying  MHD  to  a  hypersonic  vehicle  is  easy  to  see  when  we  compare  the  best 
representative  of  a  real  hypersonic  vehicle  power  system:  the  Shuttle  orbiter.  The  Shuttle  uses  fuel  cells  and 
hydrazine  powered  auxiliary  power  units  for  electrical  and  hydraulic  power.  The  three  fuel  cells  of  the  Shuttle  are 
capable  of  a  maximum  continuous  output  of  21  kilowatts  with  15-minute  peaks  of  36  kilowatts  with  a  255  lbs 
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package.  The  typical  average  power  consumption  of  the  orbiter  is  14  kilowatts  (7  kilowatts  are  available  for 
payloads).  The  hydrazine  APU’s  provide  the  hydraulic  power  for  Shuttle  Orbiter  during  the  ascent,  descent,  and 
landing.  The  APU’s  power  the  aero  control  surfaces;  engine  thrust  vector  control,  engine  valves,  landing  gear, 
brakes,  and  steering.  The  APU’s  typically  operate  for  90  minutes  over  a  power  range  of  6  kW  to  1 10  kW.  Proposals 
for  an  electrical  APU  replacement  estimate  the  units  would  weight  80  lbs.  If  fuel  cells  were  used  instead,  the  peak 
1 10  kilowatts  for  the  equivalent  replacement  electrical  actuators  in  the  orbiter  would  require  1275  lbs  of  fuel  cell  to 
provide  the  total  potential  electrical  power  in  the  entire  vehicle.  The  initial  studies  of  HVEPS  show  that  an  MHD 
generator  system  (magnet,  channel,  seeding)  can  produce  a  10  megawatt  burst  for  10  minutes  in  a  package  of  less 
than  4000  lbs  where  the  equivalent  weight  for  a  Shuttle  APU  would  be  7000  lbs  and  for  fuel  cells  would  be  121000 
lbs.  Add  to  these  weights,  auxiliary  power  units  woidd  require  heavy  amounts  of  hydrazine  and  the  fuel  cells 
require  heavy  amounts  of  oxidizer  to  operate.  Because  the  MHD  generators  extract  power  directly  from  scramjet 
flow,  they  generate  the  electrical  voltage  and  current  as  soon  as  the  seeding  is  added  to  the  flow  stream,  in  near 
instantaneous  response  times  of  milliseconds,  ideal  for  many  directed  energy  weapons. 

Any  attempt  to  use  the  extremely  high  speed  gas  flows  to  power  a  turbine  woidd  either  directly  destroy  the 
internal  turbine  blades  in  the  extreme  heat  of  the  inlet  gas  or  would  require  an  unacceptably  enormous  cooling 
system  to  prevent  damage.  MHD  generators  work  optimally  at  the  high  temperatures  of  hypersonic  flight.  Studies 
have  shown  that  hypersonic  vehicle  will  likely  require  power  levels  into  100’s  kilowatts  (possibly  even  megawatts) 
of  power  for  aero  control  surface  electrical  actuators  given  the  need  for  more  maneuverability  over  flight  durations 
of  several  hours.  MHD  generators  can  directly  power  the  electrical  actuators  or  can  be  applied  to  a 
magnetogasdynamic  (MGD)  “virtual  control  surface”  to  deflect  the  natural  plasma  surrounding  the  hypersonic 
vehicle.  The  natural  (or  microwave  induced)  plasmas  can  be  channeled  through  a  MGD  “virtual  cowl”  at  the 
scramjet  inlet  to  increase  inlet  flow.  In  addition,  startup  of  the  scramjet  may  require  a  plasma  torch  to  initiate 
combustion.  Finally,  a  high  power  generator  woidd  make  hypersonic  aircraft  strong  candidates  for  directed  energy 
weapons.  The  final  test  effort  of  this  program  is  verification  of  a  MHD  computational  fluid  dynamic  computer  code 
that  has  been  developed  by  Purdue  University  and  the  University  of  Tennessee  Space  Institute. 

II.  Technical  and  Historical  Background  of  MHD 

MHD  generators  produce  a  current  in  an  ionized 
gas  by  passing  the  gas  through  a  magnetic  field.  The 
magnetic  electric  field  (U  x  B)  consists  of  a  perpendicular 
(Faraday)  component  and  an  axial  (Hall)  component  relative 
to  the  direction  of  motion.  Electrons  moving  because  of  this 
field  first  form  surface  charges  at  of  the  MHD  channel  and 
develop  an  opposing  electric  field.  Electrodes  are  positioned 
to  use  the  Faraday  or  Hall  fields,  but  diagonally  connected 
electrodes  generally  produce  the  optimum  configuration  in 
rectangular  channels  by  taking  full  advantage  of  both  the 
Faraday  and  Hall  electric  fields  simultaneously.  The 
majority  of  MHD  research  has  been  with  diagonally 
connected  generators,  although  some  pidse  power 
generators  utilized  continuous  Faraday  electrode  designs. 

Michael  Faraday  proposed  the  initial  concepts  of 
MHD  in  the  1830’s.  The  first  demonstration  of  an  electric 
motor  was  based  on  a  conducting  wire  rotating  in  the  bath  of  liquid  mercury  and  was  actually  a  demonstration  of 
MHD.  Concepts  for  electrically  conducting  gas  MHD  generators  began  to  appear  in  patents  as  early  as  1910; 
however,  no  methods  of  ionization  were  suggested,  so  actual  devices  were  never  built.  The  first  working 
experiments  on  a  conducting  gas  MHD  generator  began  with  Karlovitz  and  Halasz  in  1938  at  Westinghouse.  This 
generator  was  a  coaxial  cylindrical  design,  but  with  limited  capability  as  the  basic  plasma  physics  was  not  well 
understood  at  that  time.  As  more  accurate  plasma  physics  models  were  applied,  more  practical  generators  were 
constructed  in  the  1960’s  and  1970’s.  MHD  aerodynamic  controls  were  design  by  Krantowitz  in  1955  with  the  use 
of  MHD  for  reactive  control  of  re-entry  vehicles  using  high  natural  thermal  ionization  produced  from  the  re-entry 
drag  at  hypersonic  speeds.  A  large  Air  Force  MHD  research  program  in  the  1960’s  and  1970’s  focused  on 
supersonic  combustor  flows  and  several  studies  used  turbine  engine  exhaust  flows. 


Figure  3  Combustion  MHD  Generator 
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The  combustion-driven  MHD 
power  generation  is  the  most  mature  and 
well  understood  type  of  MHD  technology 
application.  Combustion  driven  generators 
need  internal  oxidizers,  which  can  be  a 
disadvantage  for  any  long  duration 
aerospace  mission.  Combustion  MHD  is, 
however,  competitive  with  almost  any 
energy  storage  and  rotating  generators  over 
durations  of  1  to  1000  seconds.  MHD  pulse 
power  systems  with  power  density  on  the 
order  100’s  MW/m3  have  been  well 
demonstrated  as  practical  systems  over  the 
pass  40  years.  Russian  work  in  1980’s  was 
reported  to  have  produced  an  output  of  500 
MW  at  a  power  density  approximately  600 
MW/m3.  The  theoretical  power  densities  of 
hypersonic  MHD  generators  can  be  made 
extremely  high  (100’s  to  1000’s  MW/m3). 
No  other  non-nuclear  power  technology  can 
compete  with  this  technology  for  hypersonic 
air  and  space  vehicle  power  sources. 

Thermal  ionization  is  a  byproduct 
of  the  high  collision  rate  of  the  atoms  in  the 
high  temperature  gas.  Achieving  natural 
thermal  ionization  in  molecular  gases 
requires  extremely  high  temperature.  For 
example,  ionization  of  pure  air  requires 
temperatures  over  4,000  °K.  The  traditional 
technique  of  adding  high  electron  affinity 
alkali  metals  is  easier  in  the  high 
temperature  flow.  Additional  ionization 
using  high  electric  fields  coronas,  electron 
beams,  or  high  power  microwaves  is 
possible;  the  practical  limitation  in  any  case 
is  the  power  is  required. 


Figure  5  Hypersonic  Vehicle  Engine  Power 


Figure  4  MHD  Power  System  Weight  Estimate 


Large  vehicle  scramjet  engines  will 
need  to  produce  gigawatts  of  power  to 
accelerate  and  reach  a  constant  cruise  speed, 
and  will  need  to  produce  hundreds  of 
megawatts  of  power  to  maintain  cruise 
speeds.  MHD  generators  can  easily  extract 
large  amounts  of  electric  power  from  the 
tremendous  kinetic  and  thermal  energy 
available  in  a  hypersonic  flow. 
Demonstrated  enthalpy  extraction  of  MHD 
generators  in  the  0. 5-2.0  MW  range  have 
been  at  the  order  of  1. 0-5.0  %  and  for  larger 
units  in  the  10-500  MW  range,  enthalpy 
extraction  increases  to  the  10.0-20.0  % 
range.  For  large-scale  MHD  power  systems, 
isentropic  efficiency  approaching  80%  is 
theoretically  possible.  Since  scramjets  are 
non-rotating  machines,  non-rotating  MHD 
generators  are  a  natural  match  to  produce  efficient  power. 


0  2000  4000  6000  8000  1  -104  1.2  104 

Range  in  Miles 

200,000  lbs  (Dry),  6000  mi/hr  Vehicle 

Figure  6  Hypersonic  Vehicle  Mass  Fractions 
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The  high-speed  gas  flow  produces  extremely  high  power  density  MHD  power  systems  that  are  a  small 
percentage  of  the  total  vehicle  weight  (4000  kg  power  system/ 100,000  kg  vehicle  =  4%).  This  weight  is  achievable 
even  with  helium  cooled  superconductor  magnets;  high  temperature  YBCO  superconductor  magnets  have  the 
potential  of  using  lightweight  liquid  nitrogen  (or  liquid  hydrogen)  temperature  cryocoolers. 


ITT.  Space  Charge  Voltage  Regulation  Concept 


The  measurements  show  in  Fig.  2  show  that 
the  actual  voltage  and  current  measurements  from  a 
MFID  generator  would  require  significant  filtering  to 
produce  a  smooth  output.  In  addition,  the  voltage  is 
directly  determined  by  the  gas  flow  velocity  because 
the  electric  field  in  the  generator  is  U  X  B  where  U  is 
the  flow  velocity  and  B  the  magnetic  field.  As  U 
varies,  the  only  alternatives  to  regulate  the  output 
voltage  are  adjust  the  magnet  current  to  change  the 
magnetic  field  or  perform  full  output  power 
conditioning  with  heavy  capacitor  filtering.  A 
generally  ignored  factor  is  that  the  current  is  a 
combination  of  the  conductivity  and  space  charge  as 
which  can  be  written  as  J  =  c(E  +  U  X  B)  +  pU,  where 
J  is  the  current  density,  a  is  the  conductivity,  E  the 
electric  field,  U  the  gas  flow  velocity,  B  the  magnetic 


CONCENTRIC  ClUNDRICAl.MAGNEIOmDRODrN.UllC  GENERATOR 


Figure  7  Conceptual  Cylindrical  MHD  Generator 


field,  and  p  is  the  space  charge  density.  The  electric 
fields  from  space  charge  in  a  rectangular  channel  tend  to  produce  high  electric  fields  in  the  sharp  corners  of  the 
channel.  Because  the  magnetic  field  across  the  channel  also  creates  basically  a  transverse  electric,  the  radial  electric 
fields  from  space  charge  will  disturb  the  magnet’s  electric  field  in  certain  regions,  so  space  charge  is  neutralized  as 
much  as  possible  in  a  rectangular  MELD  generator  to  minimize  this  effect. 


The  symmetry  of  the  space  charge’s  electric  field  as  calculated  by  Gauss’s  law  indicates  that  a  symmetric 
geometry  such  as  a  cylinder  or  sphere  can  take  advantage  of  the  space  charge  electric.  This  in  fact  was  the  physical 
basis  for  the  performance  of  electron  tubes  as  show  in  Fig  4  above  where  electrons  are  the  actual  space  charge.  The 
triode  in  fact  introduces  a  control  grid  to  regulate  the  voltage  between  the  cathode  and  anode.  The  electric  field  is 
easily  calculated  for  a  cylindrical  geometry  with  an  enclosed  charge  using  Gauss’s  Law.  In  the  case  of  MHD 
channel  with  space  charge,  the  enclosed  charge  is  assumed  distributed  uniformly  in  a  long  cylindrical  geometry  to 
produce  primarily  a  radial  electric  field  (the  end  effects  of  the  field  are  assumed  to  be  negligible).  The  radial  electric 
field  of  the  space  charge  aligns  with  the  radial  electric  field  created  from  the  circular  magnetic  field. 


cathode 

anode 


Precision  Graphics 


Figure  8  Concentric  Electrodes  in  Triode 


Er(Rout) 

Figure  9  Radial  Electric  Field  in  a  Cylindrical  Channel 


5 


IV.  Theoretical  Analysis  of  Space  Charge  Effects 

A.  Generator  Output  Voltage  and  Electric  Field 

We  will  now  derive  the  equations  and  plots  that  show  the  effects  of  space  charge  on  the  output  voltage, 
electrode  electric  field,  magnetic  field,  electrical  efficiency,  generator  pressure  drop,  and  extracted  enthalpy.  The 
calculations  will  demonstrate  that  adjusting  the  space  charge  is  a  reasonable  means  to  regulate  the  output  voltage 
(electric  field)  and  reduce  the  magnetic  field  required  for  a  particular  power  level.  For  the  plot  calculations,  the 
generator  power  level  is  set  at  5  megawatts  and  the  generator  dimensions  are  Router  =  1  m,  Rinner  =  0.9  m,  L  =  1  m, 
(the  resulting  power  density  is  Pd  =  8.4  MW/m3).  Assume  the  propulsion  is  provided  by  a  hydrocarbon  fueled 
scramjet.  Kerosene  and  oxygen  give  a  reasonable  representation  of  the  electron-molecular  collision  cross-section 
(approximately  Qe  =  40  x  10'2CI  m2)  for  a  hydrocarbon,  and  if  the  combustion  temperature  is  about  Tin  =  3000  °K, 
the  electron  mobility  works  out  to  be  Ke  =  0.528  m/s/V/m.  We  will  assume  the  electron  density  is  sufficient  to  yield 
an  electron  conductivity  oe  =  -20  mhos.  The  minus  sign  is  used  in  the  electron  conductivity  to  calculate  the  correct 
signs  and  values  for  voltage  and  electric  field  because  the  space  charge  fraction  must  be  5  <  0  for  a  plasma 
dominated  by  electrons  and  negative  ions.  The  inlet  gas  pressure  is  set  at  Pjn|et  =  1  atm  and  the  true  Mach  numbers 
are  determined  using  the  temperature  dependent  speed  of  sound,  Cm:=^/y-Rg-Tin*,  where  the  Universal  gas 
constant  for  air  is  Rg(air)  =  287.05  Joule/Kg/°K,  and  the  constant  pressure  specific  heat  is  Cp(air)  =  1005  J/Kg/K 
and  the  ratio  of  specific  heat  is  y  =  1.4. 

The  complete  current  density  is  determined  by  Ohm’s  Law  with  the  Flail  current  and  space  charge  terms, 

>  ->  -»  i  >  -»\  -» 

J  =  Gij-I  E+  Ux  B  +  —  J xbI  +  p-U  (1) 


where  the  current  density,  electric  field,  magnetic  field,  and  flow  velocity  vectors  are  in  cylindrical  coordinates, 


(  Jr  2 

■» 

^Er^ 

-> 

f  o  > 

-» 

(  °  ^ 

JO 

E  = 

0 

U  = 

0 

B  = 

B0 

k  Jzj 

VE  z) 

,Uz) 

v  0  J 

The  magnetic  field  Be  between  concentric  cylinders  with  inner  radius  Rjnner  and  outer  radius  Router, 


B0(r)  = 


pO-N-I 

2-jr-r 


(3) 


The  maximum  magnetic  field  occurs  at  Rjnner 

pO-N-I 

Bm=  — - 

2-n-Rin 

The  internal  magnetic  field  can  be  defined  with  the  maximum  magnetic  field, 

Rin 

B0(r)  =  Bm - 

r 

Conductivity  tensor  for  electrons  (or  ions)  is 

a  ..  -ct-(K-BO) 


i  +  (kbo)“  i  +  (k-bo)' 


(4) 

(5) 


(6) 
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The  conductivity  is  the  product  of  mobility  and  charge  density,  <re  =  Ke-pe  and  ai  =  Ki-pi.  The  total  current  density, 


J  =  K 


oeijj-lE  +  U  x  B)  + 


™ij 


ae 


■  A 


V  Pi 


ij 

pe  J 


■(j  x  b)  +  (pi  -  pe)-U 


(7) 


The  electron  mobility  (Ke)  and  the  ion  mobility  (Ki)  from  the  electron  collision  frequency  (ve)  and  ion  collision 
frequency  (vi)  are  determined  by, 


e  1 

Ke  = - 

me  ve 


(8) 


Ki=  Z- 


mi  vi 


Where  the  electron  charge  is  e  =  1.60217733*10'^  coulombs  and  the  electron  mass  is  me  =  9. 1093897*  10”^^  kg. 
The  ion  mass  is  m;  =  A-AMU  where  A  is  the  atomic  weight  and  AMU  =  1.6605402*10"^  kg  (Atomic  Mass  Unit) 
and  the  ion  charge  number  (or  the  electrons  removed)  is  Z.  The  electron  collision  frequency  (ve)  and  the  ion 
collision  frequency  (vi)  are  related  to  the  electron-molecular  and  ion-molecular  collision  cross-sections  (Qe,  Q;), 


ve  =  nO-Qe- 


vi  =  nO-Qi- 


l-k-  Tin 


71-me 


16  k- Tin 


7t- mi 


(9) 


where  the  gas  density  is  n0  =  Piniet/k-Tiniet  or  the  ideal  gas  law  at  inlet  temperature  and  gas  pressure  (Boltzman 
constant  k  =  1.380658*  10“^^  J/K).  The  electron  mobility. 


Ke  = 


me 


nO-Qe- 


8-k-Tin  me  Pin 


71-me 


k-Tin 


•Qe- 


8-k-Tin 


71-me 


(10) 


The  ion  mobility 


Ki  =  Z- 


mi 


n0-Qi- 


16-k-Tin 


7i- mi 


mi  pin 


k-Tin 


Qi- 


16-k-Tin 


7t- mi 


(11) 


The  ratio  of  electron  mobility  to  ion  mobility  is  60.38(Qi/Qe)A‘ 2  and  since  Qi  »  Qe,  the  electron  mobility  is 
substantially  larger  than  any  ion  mobility  (Ke  »  Ki).  The  greater  mobility  of  the  electron  allows  the  current  density 
to  be  accurately  approximated  (in  matrix  form), 


1 


A  Jr  A 

J0 

(  iz) 


1  +  (Ke-Be)“ 
0 

Ke-B0 


0 


-Ke-B0 


1  +  (Ke-Be)' 
1  0 
1 


0 


1  +  (Ke-B0)“  1  +  (Ke-Be)’ 

The  space  charge  term  is 


A  ErA 
0 

VEV 


(  n  A 


(  n  \ 


B0 


-Ke- 


A  Jr  A 

J0 

(  Jzy* 


A  0  A 

B0 

V  0  J. 


or  ae 
Ki  Ke 


A  0  A 

0 


(12) 


pi-  pe  : 


ai 

Ki 


oe 

Ke 


(13) 


The  neutral  charged  plasma  condition  is  pi  -  pe  =  0,  and  the  ionized  plasma  conditions  are  pe  >  pi  (net  negative 
space  charge)  and  pe  <  pi  (net  positive  space  charge).  The  space  charge  densities  can  be  expressed  as  a  fraction  of 
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the  gas  density  nO,  a-e-nO  =  pe  (positive  ion  fraction  -  a),  P-e-nO  =  pi  (electron  fraction  -  P),  x-enO  =  Pn  (negative 
ion  fraction  -  /).  The  net  space  charge  becomes, 


pi  -  pe  =  a-e-nO  -  p-e-nO  =  (a  -  p)-e-nO  =  — — — -pe  =  &—  (14) 

P  Ke 

where  5  =  (a  -  P)/a  is  the  relative  excess  fraction  of  electron  charge.  The  plasma  created  by  thermal  ionization  is  a 
loosely  bound  cloud  of  electrons  and  ions  which  appear  from  the  outside  to  be  neutral.  A  net  space  charge  of 
electrons  can  be  added  externally  to  the  plasma  by  an  electron  beam  and  a  net  space  charge  of  negative  ions  can  be 
added  by  a  negative  corona.  A  complete  mixture  of  free  electrons  (a),  positive  ions  (P),  and  negative  ions  ('/_)  can 
be  described  by  a  similar  term,  A, 


pi  -  pn  -  pe 


/  \  a  -  y  -  P  ae 

=  (a  -  y  -  p)-e-nO  = - pe  =  A 


P 


Ke 


(15) 


where  A  =  (a  +  %  -  P)/a  is  the  excess  fraction  of  electron  charge  (with  negative  and  positive  ion).  The  possible 
charge  conditions  are:  neutral  space  charge  5  =  0  or  A  =  0;  negative  space  charge  5  >  0  or  A  >  0;  positive  space 
charge  5  <  0  or  A  <  0. 

The  electron  dominated  approximation  for  total  current  density  using  the  excess  charge  term  (5), 


f  Jr> 
J0 

y  Jzy 


1 


1  +  (Ke-B0)“ 
0 

Ke-B0 


-Ke-BQ 

1  +  (Ke-B@)2 
0 
1 


1  +  (Ke-B0)“  1  +  (KeB@)~ 

In  matrix  form,  the  current  density  is, 


^Er^ 

f  °  ^ 

f  0  ] 

^  Jr^ 

r°] 

(  0] 

-ae- 

0 

+ 

0 

X 

B0 

-Ke- 

J0 

X 

B0 

+  f  5-— Y 

0 

V  Ke] 

vEzy 

vUzy 

v  0  y 

l  Jzj 

v  0  j 

yUzy 

(16) 


f  Jr  ^ 
J0 

V  Jz  J 


"[(l  +  §)-Uz-B0  -  Er]  •ae 
0 


(Uz-8  -  Ke-Ez) 


ere 

Ke 


(17) 


The  three  load  connections  possible  are  the  transverse  Faraday  load  (Ez  =  0),  the  axial  Hall  load  (Er  =  0)  and  the 
diagonal  load  (non  zero  Ez  and  Er).  For  a  Faraday  connection,  the  axial  electric  field  is  Ez  =  0, 

"[(l  +  5)-Uz-B0  -  Er]  •ae 
0 


f  Jr\ 
J0 

Wz] 


(Uz-8) 


ae 

Ke 


(18) 


The  power  density  is  determined  using  the  total  electric  field  applied  to  the  electrons, 

"[(l  +  §)-Uz-B0  -  Er]  •ae 
0 

(iJz-  §)  •  — 

Ke 


Pd  =  J-(E  +  U  x  B) 


r  Er^ 

f  0  ^ 

f  °  ^ 

0 

+ 

0 

X 

B0 

U  ; 

luzy 

l  o  j 

(19) 


Performing  the  multiplication 

Pd  =  -ae-(lJz-B0  -  Er)-[(l  +  5)-(Uz-B0)  -  Er] 


(20) 


This  implies  the  minimum  conductivity  is, 


(b0-Uz  -  Er)-[(l  +  §) •  (iJz- B0)  -  Er] 
The  Hall  load  generator  current  density  becomes 

,TN  H(l  +  s)-Uz-B0]-oel 


JG  = 


Jz  J  (Uz-8  -  Ke-Ez)- 


The  Hall  generator  electrical  power  density, 
~[(l  +  5)-UzBG]  •ere 
0 

Pd  = 

(Uz-8  -  Ke-Ez)--^- 


0  I  I  0  )  I  0 

0  +  0  x  BG 

EzJ  I  Uz  J  0 


Performing  the  multiplication, 


Pd  =  -ere 


(l  +  S)-(Uz-BG)7  +  |l  -8- 


Ke-Ez 


The  diagonal  generator  current  density  (Er  and  Ez  are  non-zero), 

f[(l  +  s)-Uz-BG-Er]-ael  r,  , 


Pd  =  JE: 


(Uz-8-  Ke-Ez) •- 


Er  1  f  0  1  f  0 
0  +  0  x  BG 

EzJ  V  Uz  J  0 


Performing  the  multiplication 

Pd  =  -ae-  1  -  (2  +  §)■  UZ'B6  -Er2  +  (l  +  s)-(Uz-BG)2  +  f  1  -  5--^-l-Ez2  (26) 

_[  Er  J  y  Ke-Ez  J 

The  total  radial  electric  field  Er  is  the  sum  of  the  fields  produced  by  the  surface  charge  of  the  electrodes,  the 
space  charge,  and  the  magnetically  induced  electric  field 

E(r)  =  Es(r)  +  Ep(r)  +  EB(r)  (27) 

The  induced  electric  field  from  the  magnetic  field  B0(r)  =  Bm  (Rinner/r) 

Rin 

EB(r)  =  Uz-BG(r)  =  Uz-Bm - 

r  (28) 
The  magnet’s  electrical  potential  is  determined  from  the  its  electric  field 


— VB  =  EB(r) 
dr 

Integrating  from  the  inner  radius  to  some  point  inside 


VB(r)=  Uz-Bm- - dR 

R 
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Using  the  change  of  variable  X  =  r/Rinner,  r  =  RinnerX,  dr  =  RinncrdX, 


VB(r)  = 


•Xr 


Rin 

Uz-Bm- - dX 

X 


The  potential  from  the  magnetic  field 

VB(r)  =  Uz-BmRin-ln 
Integrating  from  the  inner  radius  to  a  point  inside 
VB(Rout)  =  Uz-BmRin-ln 


Rout  ) 

Rin  J 


(31) 


(32) 


(33) 


The  electric  field  from  the  electrode’s  surface  charge  in  a  classic  rectangular  MHD  channel  is  a  constant 
value  and  so  the  potential  is  a  simple  function  Vs  =  x  Es,  where  Es  is  the  surface  charge  electric  field  and  x  is  the 
width  of  the  channel.  Since  the  surface  charge  is  induced  by  the  magnet  electric  field,  the  electrode  electric  field  is 
usually  defined  by  a  proportional  constant  K  called  the  “load  factor”  so  that  Es  =  K  Eb.  This  is  unique  to  the 
rectangular  geometry  since  load  factor  is  correctly  defined  by  the  surface  charge  potential  and  the  magnet’s  electric 
potential  (Vs  =  K  Vb).  The  “load  factor”  is  a  more  complicated  concept  in  a  cylindrical  geometry  since  any  voltage 
is  dependent  on  geometry,  so  that  load  factor  must  be  defined  for  at  a  specific  location. 

We  now  calculate  the  space  charge  potential  from  Gauss’s  law.  The  electrode  surface  charge  density  2 
over  an  electrode  surface  area  (radius  Rin  and  length  Lelec)  acts  like  a  space  charge  ps 


2-7t-Rin-Lelec 

ps  =  2 - 

2 

jtRin-- Lelec 


Rin 


(34) 


The  potential  from  the  surface  charge  is  determined  from  Poisson’s  equation  (or  equivalently  Gauss’s  law) 


Transposing  rdr. 


1  d  (  d  1  2 

— —  r-— Vs(r)  =  2 - 

r  drl  dr  )  Rin-sO 


d|  r--Vs(r)  I  =  2-— •  — -dr 


dr 


sO  Rin 


(35) 


(36) 


Integrating  from  r  to  Rin, 


r-— Vs(r)  = 
dr 


2  )  R 


2 — 


■  dR : 


r2  2 


sO  J  Rin  Rin  eO 


- Rin — 


eO 


(37) 


Rin 


Dividing  by  r. 


dr 


r  Rin)  2 


-Vs(r)  =  — - •- 


Rin  r  )  eO 


(38) 


The  surface  charge  potential  becomes 

r  r 

Vs(r)  = 


R  Rin )  2 

- - -  — dR 

Rin  R  J  sO 


(39) 


Rin 
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Applying  the  change  of  variable  X  =  r/Rinner,  r  =  RinnerX,  dr  =  RinnerdX, 


Vs(r)  = 

The  integral  becomes 


•Xr 


Vs(r)  =  (xr  -  ln(Xr)  -  l  )-  — • 

2  sO 


X - -  —  RindX 

X  )  eO 


X  Rin 


The  result  is  the  surface  charge  voltage  inside  the  channel. 


Vs(r)  = 


4-W, 


Rin  ) 


X  Rin 
2  eO 


V  Rin/ 

The  plasma’s  space  charge  potential  is  determined  from  Poisson’s  equation  (Gauss’  Law) 

7  dr[_r  U1  +  dzdz 
The  axial  electric  field  Ez  is  determined  by 


d_d_v  _  pi-pe  §  ere 


sO 


eO  Ke 


Ez=  —V 
dz 

But  for  the  Faraday  generator  Ez  =  0 

— — V  =  —  Ez  =  0 
dzdz  dz 

The  space  charge  potential  will  be  calculated  for  a  uniform  space  charge  to  maintain  simplicity, 

d 


J_  d 

r  dr 


-Vp(r) 


dr 


5  ae 
)\  sO  Ke 


Transposing  rdr 


r-  — Vp(r) 
dr 


§  ere 

= - rdr 

eO  Ke 


Integrating  the  equation  over  r  to  Rin 


id,,  8  a?  r“  -  Rin“ 

r-  —  Vp(r)  = - 

dr  )  sO  Ke  2 


Dividing  by  r 


The  integral  can  be  rewritten  as, 
rRout 


(d  > 

1  -  5 

ae 

(  2\ 

Rin 

— Vp(r) 

r - 

\dr  j 

!  2e0 

Ke 

^  r  ) 

/-Rout 


dr 


Vp(r)  dr  = 


Rin 


5  as  Rin 
eO  Ke  2  1  Rin 


Rin 

r 


dr 


Rin 


(40) 


(41) 


(42) 


(43) 


(44) 


(45) 


(46) 


(47) 


(48) 


(49) 


(50) 
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Again  we  apply  the  change  of  variable  X  =  r/Rinner  (Xr  =  l/Rin),  r  =  RinnerX,  dr  =  RIIII1L.rdX, 


Vp(r)  = 


•Xr 


X 


O  8  oe  Rin“  if  8  cre^)  .  2  fV  2  „ 

— - dX  =  —  ■ - -Rm  -^Xr  -  2- 


X)  eO  Ke  2 


4  {  eO  Ke 


•  ln(  Xr)  -  1 


(51) 


The  space  charge  potential  becomes 


Vp(r)  = 


1  I  8  cse  1  2 

Rm 


4  1  eO  Ke 


r 

Rin  J 


Y  ( 

-2- In  - 


V  Rin 


-  1 


The  total  potential  is  the  sum  of  the  surface  charge,  space  charge,  and  magnetic  field  potentials 
V(r)  =  Vs(r)  +  Vp(r)  +  VB(r) 

Substituting  the  various  potentials,  the  complete  potential  is, 


(52) 


(53) 


V(r)  = 

The  total  load  voltage 


X  Rin )  1  (  8  ere  )  .2 

+  - I -Rin 

2  sO 


4  1  eO  Ke 


_l)  _2,nfj_)_, 


Vload  = 


X  Rin 
2  eO 


1(8  ae  )  .2 

Rin 


4  l  eO  Ke 


f  RoutV 


-2- In 


Rin  J 


Rout 


VRiny 


+  Uz-BmRin-ln 


-  1 


+  Uz-BmRin-ln 


f  Rout 


V  Rin 


V  Rin  J  y  Rin 

From  the  total  load  voltage  the  induced  surface  charge  on  the  inner  electrode  can  be  calculated 

Rout 


(is) ,54) 


(55) 


X  =  2- 


eO 

Rin 


Vload  -  Uz-BmRin-ln 


Rin 


RoutV 
Rin  J 


2- In 


Rout 

Rin 


1  I  8  oe  ]  .2 

- - l-Rin 

4  eO  Ke 


1 


(56) 


The  traditional  “load  factor”  definition  is  K  =  Vload/VB(Rout) 

X  Rin 
2  eO 


K: 


1  (  8  ere  |  2 

+  - - l-Rin 

4  l  eO  Ke 


Uz-BmRin-ln 


f  Rout 


l  Rin 


f  Rout 


.V  Rin 


«)  -2.i„VeV 

in  J  \  Rin  J 


+  1 


(57) 


Load  factor  is  now  dependent  on  the  space  charge  fraction  as  well  as  the  ratio  of  Rout/Rin.  The 
dependence  on  space  charge  factor  means  the  load  factor  changes  as  space  charge  changes.  We  can  apply  a 
condition  that  if  the  voltage  contributed  by  the  surface  charge  is  negligible  compared  to  the  voltage  from  the  space 
charge,  (Vs  «  Vp  which  is  the  same  as  X  =  0),  then  the  internal  voltage  of  the  generator  strictly  becomes  a  function 
of  the  magnet’s  electric  field  and  the  space  charge  of  the  plasma  (the  outer  electrode  surface  charge  does  not 
contribute  to  the  internal  electric  field  because  of  Gauss’s  Law).  The  condition  that  space  charge  voltage  is  greater 
than  the  surface  charge  voltage  is, 


1  > 


X  Rin 
sO  2 


"l  1 

r  § 

ere  ^  .2 

— 

- -Rin 

4  1 

ysO 

Ke  ) 

(58) 
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The  lower  limit  of  space  charge  fraction  to  maintain  space  charge  voltage  greater  than  surface  charge  voltage  is, 


Vload  -  UzBmRinl 


5>  - 


V  Rin  ) 


f  RoutV 

.1  Rin  J 


-  2- In 


(  RouV 


V  Rin  ) 


-  1 


(59) 


-■[  l-Rin2 
2  l  eO  Ke 


This  becomes  the  dividing  line  between  classic  magnetohydrodynamics  and  electromagnetohydrodynamics,  where 
the  independent  space  charge  electric  field  dominates  the  electrode  surface  charge  electric  field  induced  by  the 
magnet’s  U  X  B  electric  field.  Assuming  the  surface  charge  is  negligible,  2  =  0, 


Vtotal(r)  = 


The  load  voltage  at  Rout, 


1  [  5  oe 


4  l  s0  Ke 


■Rin2' 


Rin  J 


)2  -  2*ln(— — 


^  Rin 


+  Uz-BmRin  lnf  - 

y  Rin 


vio,d=i  2  -  W 

4  l  eO  Ke 


f  RoutV 


•Inf 


i  Rm  J 

The  generator  voltage  is  controlled  by  the  space  charge  fraction. 


Rout 

Rin 


+  Uz-Bm  Rin- 


di-  In 


Rout 


V  Rin 


(60) 


(61) 


Space  Charge  Fraction 

Figure  10  Generator  Output  Voltage 
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The  maximum  magnetic  field  can  be  determined  from  the  load  voltage  equation. 


Space  Charge  Fraction 

Figure  11  Magnetic  Field  can  be  reduced  as  Space  Charge  changes 


The  magnet  current  is  determined  using  the  definition  of  the  maximum  magnetic  field 


Bm  = 


pO-N-I 

2-jt-Rin 


Solving  for  current  and  substituting  the  value  for  the  maximum  magnetic  field,  the  magnet  current  becomes 


2n  Rin 
pO  N 


1  +  2- In 


Rout 

Rin 


Rout^j 

Rin  J 


2 


oe  ] 

- Rin 

eO-Ke  J 


Vload 

H - 

Rin 


Uz-  In 


Rout^j 

Rin  J 


(63) 


(64) 
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The  magnet  current  needed  decreases  as  space  charge  fraction  increases, 


Space  Charge  Fraction 

Figure  12  Magnet  Current  per  Magnet  varies  with  Space  Charge 

The  space  charge  fraction  can  also  be  determined  from  the  load  voltage  equation. 


5  =  4-sO- 


Ke 


Uz-Bmln 


Rout 
\  Rin 


Vload 

Rin 


ere 


Rin 


1  +  2- In 


f  RouA  f  RoutV 


V.  Rin  J  V  Rin  J  _ 

The  output  voltage  and  magnetic  field  can  be  regulated  by  space  charge  fraction. 


Voltage  Regulation 


Gas  Flow  Velocity  (m/s) 

Figure  13  Voltage  Regulation  using  Space  Charge 


(65) 
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The  space  charge  current  is  determined  by  the  equation  Is  =  8*(ce/Ke)*Uz*ji(Rout  -  R i n  ) 


0  1000  2000  3000  4000  5000  6000  7000  8000  9000  1  -104 


Gas  Flow  Velocity  (m/s) 


Figure  14  Space  Charge  Current  regulates  Voltage 


The  radial  electric  field  in  the  channel  with  non-zero  electrode  surface  charge 


If  r  Rin)  Z  1  5  ere  f  r  Rin  )  Rin 

Etotal(r)  =  — •  2 - - 1 - Rin- - +  Uz-Bm - 

2  y  Rin  r  )  s0  2  e0  Ke  V  Rin  r  )  r 

The  total  electric  field  with  the  inner  electrode  grounded, 

1  5  oe  (  r  Rin^i  Rin 

Etotal(r)  = - Rin- - +  Uz-Bm - 

2  s0  Ke  l  Rin  r  )  r 


(66) 


(67) 


The  electric  field  at  the  outer  electrode  as  space  charge  fraction  changes, 


Space  Charge  Fraction 

Figure  15  Generator  Electric  Field  at  Outer  Electrode  Surface 
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To  determine  the  relationship  for  5  with  a  constant  power  density.  Starting  from  the  Faraday  power  density 


Pd  =  -(Uz-B0  -  Er)-[(l  +  s)-Uz-B0  -  Er]-oe 


(68) 


We  will  define  the  term  Ed  so  that  Er  =  5  Ed 


Pd  =  -(Uz-B0  -  S-Ed)-[(l  +  S) •  Uz- B0  -  S-Ed]-cre 


(69) 


The  result  is, 


Pd  = 


(uz-B0-Ed  -  Ed2)-82  +  (-Uz2-B02  +  2-Uz-B0-Ed)- 


?  2 
•S-Uz“B0 


oe 


The  terms  in  parentheses  can  be  defined  as 


.  =  ( 


Pd  =  \Al-8  +  BI  S  -  Cl/cse 

A1  =  Ed- (Uz-B0  -  E8)  . B1  =  -Uz-B0-(Uz-B0  -  2- Ed)  ;C1  =  (Uz-BO)' 
the  solutions  for  space  charge  fractions  in  Equ.  67  are, 
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— R1  1  f  R1  1  Prl  r  1 
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f+  1 

Pd 
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4 
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J  A1 

ae 
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2 

.  'I 

^  ~B1 
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.2 
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^  B1 

J  A1 

Pd  Cl 

oe  Al_ 

4 

vaT, 

The  space  charge  fraction  for  constant  output  power, 


(70) 


(71) 


(72) 
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Figure  16  Power  Regulation  with  Space  Charge  Changes 
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The  space  charge  current  for  constant  output  power 
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Figure  17  Power  Regulation  as  Space  Charge  Current  changes  with  Gas  Flow  Velocity 


B.  Generator  Efficiency 


The  local  efficiency  pe  is  defined  by  the  ratio  of  electrical  power  density  to  mechanical  power  density 


rje  = 


I  Pd  I 

|-»  >\ 

|u-f| 


(73) 


Since  we  know  the  electrical  power  density,  we  now  need  to  determine  the  mechanical  power  density  U*f  from  the 
force  density, 

"[(l  +  5)-Uz-B0  -Er]-cre~ 

0 


>  ->  >  -> 
f=  p-E+  Jx  B=  §■ 


^Er^ 

„  oe 

§ - 

0 

Ke 

s0y 

(Uz-S)-— — 
Ke 


f  0  A 


B0 

voy 


(74) 


This  simplifies  to, 


(Er-Uz-B0)-§-  — 

Ke 

0 

[(§+  l)-Uz-B0  -  Er]-B0  •aej 


(75) 


The  mechanical  power  expended  on  the  conducting  gas  flowing  with  speed  U, 


Uf=  U-vp -E  +  Jx  Bj 


Substituting  the  vectors  for  the  velocity,  electric  field  and  magnetic  field 


->  > 
U-f  = 


(0A 

0 

'vUz; 


(Er-Uz-B0)-§- 


ae 

Ke 


[(§+  l)-Uz-B0  -Er]-B0-ae_ 


(76) 


(77) 
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Performing  the  multiplication 

U-f  =  [(l  +  S)-Uz-B0  -  Erj-Uz-B0-cse 
Now  we  can  calculate  the  electrical  efficiency 
Pd 
U-f 


(78) 


(79) 


Substituting  the  electrical  power  density  and  mechanical  power  density 
_  [-(Uz-B6  -  Er)-[(l  +  s)-Uz-B0  -  Erjj-cre 
116  [(l  +  5)-Uz-B0  -  Er)-Uz-B0-cre 

The  electrical  efficiency  simplifies  to, 


(80) 


-Uz-B0  +  Er 
pe  = - 

Uz-B0  (81) 

Applying  the  space  charge  electric  field  and  the  magnetic  induced  electric  field,  the  electrical  efficiency  become, 

1  8  ere  Rout  f  Rout  Rin  ^ 

rie= - - -1  (82) 

2  s0  Ke  Uz-Bm  y  Rin  Rout ) 


Space  charge  fractions  for  a  constant  efficiency 

pe  +  1 


Ke  Uz-Bm 

8  =  2e0 - 

ere  Rin 


f  Rout  V 

.1  Rm  J 


-  1 


The  space  charge  fraction  for  a  constant  generator  electrical  efficiency  is  shown  in  Fig  17. 


(83) 
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Figure  18  Electrical  Efficiency  Regulation  as  Space  Charge  Fraction 
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The  space  charge  current  for  a  constant  generator  electrical  efficiency 
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Figure  19  Electrical  Efficiency  Regulation  with  Space  Charge  Current 

Local  efficiency  for  constant  load  voltage  and  magnetic  field 
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^  Rin 


Rin 
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1  +  2- In  - 


f  RouA  f  RoutV 


^  Rin 


l  Rin  )  _ 


(84) 


The  local  or  electrical  efficiency  is  the  only  efficiency  directly  affected  by  changes  in  the  electromagnetic  fields.  All 
other  efficiencies  are  primarily  affected  by  thermodynamic  properties  and  only  indirectly  by  the  electrical  properties 
through  the  dependence  on  the  local  efficiency. 

The  isentropic  efficiency  depends  on  the  local  efficiency, 

pe 


pp  = 


1  + 


T-O-pe) 
2  A  Cm  J 


y  - 1 


(85) 


The  total  generator  efficiency  in  turn  is  calculated  from  the  isentropic  efficiency, 

J1P  •(•/->) 

l-(l-Rp) 


y 


fig  = 


y-i 

y 


(86) 


l-(l-Rp) 

The  generator  efficiency  depends  the  thermodynamic  cycle  properties  of  the  generator,  specifically  the  output  to 
input  pressure  ratio  (Rp)  and  the  specific  heat  ratio  (y  =  Cp/Cv). 

The  space  charge  dependence  of  the  generator  efficiency  is  determined  by  solving  Eq  76  for  pp, 


-ln(Rp) 


Y-l 


-ln(Rp) 


y-l 


pp  = 


Iny-pg-e 


+  pg  +  e 


■y  +  ln(Rp)-y  -  ln(Rp) 


(87) 


ln(Rp)-(y  -  l) 
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In  turn  we  solve  Eq.  75  for  qe, 


This  is  substituted  into 


The  plot  of  space  charge  fraction  for  a  constant  generator  electrical  efficiency 


(88) 


(89) 
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Figure  20  Generator  Efficiency  with  Space  Charge  Fraction 

The  plot  of  space  charge  current  for  a  constant  generator  electrical  efficiency, 
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Figure  21  Generator  Efficiency  dependency  on  Space  Charge  Current 
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C.  Internal  Pressure  Drop  and  Enthalpy  Extraction  of  Generator 

The  axial  force  density  for  the  Faraday  generator  is, 

fz=  [(5  +  l)-Uz-B9  -  Erj-B9-cre 
Using  the  minimum  conductivity  for  the  Faraday  generator 

-Pd 

ere  = 


(b9-Uz  -  Er)-[(l  +  §)-(Uz-B9)  -  Er] 
The  axial  force  density  becomes 

fz=  [(§+  l)-Uz-B9  -  Er]-B9- 


-Pd 


(b9-Uz  -  Er)-[(l  +  §)-(Uz-B9)  -  Er] 


Simplifying  the  equation 
fz  =  Pd- 


B9 


Er  -  Uz-B9 

The  axial  pressure  drop  is  calculated  from  the  integral  of  the  force  density  over  the  length  of  the  generator. 


APz  : 


■L 


Pd- 


B9 


dz 


0 


The  integral  becomes 

APz=  Pd- 


Er  -  Uz-  B9 


B9 


Er  -  Uz-B9 
Using  the  value  for  the  electrical  efficiency 

Er  -  Uz-  B9 


Uz-qe  = 


B9 


The  axial  pressure  drop  becomes 

„  Pd  L 

APz  - - 

qe  Uz 

The  electrical  efficiency  can  be  written  in  terms  of  the  constant  pressure  drop 


*  „  A  Ke  Uz-Bm 

8  =  2e0 - 
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Uz  APz 


f^|  -1 
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(90) 
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The  plot  of  generator  pressure  drop  with  space  charge  fraction 
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Figure  22  Pressure  Drop  Regulation  with  Space  Charge 

The  plot  of  space  charge  current  for  constant  generator  pressure  drop, 
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Figure  23  Pressure  Drop  Regulation  with  Space  Charge  Current 


The  reduction  of  thrust  caused  by  the  generator  is 
Th  =  Area-APz 


Space  charge  fraction  for  constant  thrust  reduction 


17  „  „  Ke  Uz-Bm 

5  =  2e0 - 

cse  Rin 
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Uz  Th 


( Rout  y 

A  Rm  J 


(99) 
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The  plot  of  constant  reduction  of  thrust  from  the  generator  is  shown  in  Fig.  24. 


Thrust  Regulation 


Figure  24  Thrust  Reduction  Regulation 

The  space  charge  current  for  a  constant  generator  thrust  reduction. 
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Figure  25  Thrust  Reduction  Regulation  with  Space  Charge  Current 


The  generator  enthalpy  drop  is  determined  by  the  equation 


Tin  Pd  L 

h  =  Rgas - - 

Pin  1 


qe 


Uz 


(101) 


The  main  byproducts  of  scramjet  combustion  have  essentially  the  same  thermodynamic  properties  of  air,  so  using 
the  universal  gas  constant  for  air  (Rgas(air)  =  287.05  Joule/Kg/K)  as  well  as  the  heat  capacity  for  air  (Cp(air)  =  1005 
J/Kg/K)  are  accurate  enough  to  calculation  of  extracted  enthalpy. 
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The  electrical  efficiency  for  a  constant  enthalpy  from  Equ., 


pe  =  Rgas-Tin-Pd- 


h-Pin-Uz 

Space  charge  fraction  for  a  constant  enthalpy  becomes, 


^  „  n  Ke  Uz-Brn 

o  =  2e0 - 

cse  Rin 


Rgas-Tin-Pd- 


h-Pin-Uz 


+  1 


Rout  V 


Rin 


J-1 


The  constant  enthalpy  is  shown  in  Fig  26. 
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Figure  26  Enthalpy  Regulation 

The  space  charge  current  for  a  constant  generator  enthalpy  is  shown  here. 
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Figure  27  Enthalpy  Regulation  with  Space  Charge  Current 
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D.  Electromagnetic  Reynolds  Number 

The  Reynolds’  number  is  generally  interpreted  as  the  ability  of  the  fluid  dynamics  to  modify  the 
electromagnetic  field.  The  definition  of  the  electromagnetic  Reynolds’  number  is  the  ratio  of  extracted  electrical 
energy  to  total  stored  energy  in  the  electromagnetic  field. 


R=  — 
ub 

The  electrical  energy  density  extracted  from  the  gas  flow  is, 

_  J_  P  d 

pe  Uz 

The  energy  density  stored  in  the  electromagnetic  field  is, 

1  B0(Rout)2  1  2  1  2 

ub  = - t - sO  Eif  Rout)  H - eOEzfRout) 

2  po  2  2 

Rewriting  the  last  two  terms  of  the  stored  energy, 

2 

,  1  B0  1  eO-  po  2  1  sO'iio  2 

ub  = - H - Eif  Rout  )  + - E^Rout) 

2  po  2  po  2  po 

Permittivity  and  permeability  are  related  to  the  speed  of  light  (e0  po  =  1/c2) 
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ub  _  l  BO  (Rout)  |  1  1 


po 


2  1  1 

ErfRout)  + - 

2  2  2 


EzfRout)" 


po-c  po-c 

Substituting  the  values  for  the  magnetic  field  and  the  Faraday  electric  field  condition  (Ez  =  0) 

,  2 

Km 

Bm- -  I  ,  .2 


ub  = 


1 


Rout  J  1  1  f  Erf  Rout)  ) 
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Now  substituting  the  value  for  the  radial  electric  field 
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The  stored  energy  density  term  can  be  written  as, 
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Since  the  velocity  ratio  Uz/c  «  1  and  the  electromagnetic  Reynolds’  number  becomes, 
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The  electrical  efficiency  for  a  constant  Reynolds’  number 
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Space  charge  fraction  for  constant  Reynolds's  number  becomes 

.2  pO 
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The  space  charge  fraction  for  constant  Reynold’s  number. 
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Figure  28  Reynolds’  Number  Regulation 

The  space  charge  current  for  constant  Reynold’s  number, 
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Figure  29  Reynolds’  Number  Regulation  with  Space  Charge  Current 
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E.  System  Weight  Benefit  from  Eliminating  the  System  Power  Conditioner 

A  system  weight  benefit  comparison  can  be  made 
between  using  space  charge  regulation  power  supply  and  a 
full  output  power  conditioner.  The  power  conditioner  is 
required  when  the  generator  voltage  does  not  match  the 
various  load  voltage.  Power  conditioners  are  notorious  high 
mass  items.  The  biggest  limitation  is  the  power  density  of 
converter  which  is  typically  on  the  order  of  lkW/kg. 

However,  for  our  analysis,  we  will  assume  that  a  large 
technology  leap  occurs  in  power  conditioners  and  that  the 
power  density  is  5kW/kg  so  that  the  power  conditioner  mass 
is  Mpc  =  1000  kg  for  a  5  MW  power  supply  (P  =  5  MW). 

The  power  conditioner  efficiency  is  assumed  to  be  qpc  = 

98%  with  the  same  large  technology  leaps.  Assume  the 
vehicle  is  moving  at  velocity  at  about  Mach  6  (Uvehicie  = 

2000  m/s)  and  the  power  system  operates  continuously  over  Figure  30  Diffuse  Electron  Beam 

the  total  mission  time  of  4  hours  (T  =  14400  s).  This  duration  is  reasonable  if  electrical  power  is  required  for  MGD 
engine  compression,  plasma  engine  ignition,  MGD  surface  control,  etc. 

The  space  charge  can  be  provided  by  a  diffuse  electron  beam  to  inject  free  electrons  or  a  high  voltage 
corona  to  inject  negative  ions  into  the  gas  flow.  Space  charging  is  not  an  attempt  to  increase  conductivity  by 
detaching  electrons  from  the  gas  atoms.  Space  charging  injects  free  electrons  or  attaching  electrons  to  a  small 
number  (space  charge  fraction)  of  gas  atoms.  The  energy  cost  of  attaching  electrons  measured  by  the  atom’s 
electron  affinity  is  in  the  range  of  1-3  eV.  Increasing  electron  conductivity  through  non-equilibrium  ionization 
requires  fully  detaching  electrons  measured  by  the  atom’s  ionization  potential  which  ranges  from  6-300  eV. 

The  power  conditioner  total  energy  cost  is, 


j  ^ 

TEpc  =  —-Mpc -Uvehicie-  +  (l  -  qpc)-P-T 


(115) 


The  total  energy  cost  to  fly  a  full  output  power  conditioner, 
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Figure  31  Total  Fuel  Energy  Required  to  fly  a  Full  Output  Power  Conditioner 
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The  fuel  mass  required  to  fly  a  hill  output  power  conditioner 
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Figure  32  Fuel  mass  Required  to  fly  a  Full  Output  Power  Conditioner 


The  minimum  current  for  space  charging  is, 
cje-Uz-L 


Ieb  =  §•- 


Ke 


(116) 


If  the  electron  beam  voltage  is  Veb  and  the  electron  beam  array  efficiency  is  qeb,  the  supply  power  is, 
Veb-Ieb 


Peb 


qeb 


The  total  energy  cost  of  the  space  charger  power  supply  is, 
1  2 

TEeb  =  — -Meb-Uvehicle  +  Peb-T 
2 

The  total  energy  cost  for  a  space  charger  power  supply, 
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Figure  33  Total  Energy  Required  for  a  Space  Charger  Power  Supply 
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The  fuel  mass  required  to  fly  the  space  charger  power  supply, 
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Figure  34  Fuel  Mass  Required  for  a  Space  Charger  Power  Supply 


The  system  weight  benefit  is  the  ratio  of  space  charger  power  supply  energy  cost  to  power  conditioner  energy  cost, 

TEeb 

TEpc  (H9) 


SystemBenefit : 


The  system  benefit  ratio  for  various  generator  output  voltages, 
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Figure  35  System  Benefit  Ratio 

The  system  benefit  curves  show  that  the  mass  of  a  space  charger  power  supply  is  only  a  small  fraction  of  a 
full  output  power  conditioner.  The  energy  cost  of  lifting  the  mass  of  a  power  conditioner  can  be  also  measured  by 
the  mass  of  hydrocarbon  fuel  (assuming  the  energy  density  for  hydrocarbons  is  about  40  MJ/kg)  needed  to  lift  a 
power  conditioner.  A  5MW  power  conditioner,  based  on  a  theoretical  5kW/kg  power  electronic  power  density, 
weights  1000  kg  or  about  the  same  as  the  projected  1000  kg  for  a  5MW  MHD  generator.  The  full  output  power 
conditioner  would  require  an  additional  435  kg  (959  lbs)  of  fuel  to  reach  the  2000  m/s  (4500  mi/hr  or  Mach  6)  cruise 
velocity. 
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y. 


Conclusions 


New  methods  for  voltage  regulation  will  be  needed  to  make  lightweight  aerospace  MHD  power  systems. 
This  paper  showed  that  space  charge  can  be  a  viable  method  of  voltage  regulation  for  a  MHD  generator.  Space 
charge  would  negatively  alter  the  electric  field  in  rectangular  geometries;  however,  this  paper  shows  that  space 
charge  has  advantages  in  cylindrical  generator  geometries.  This  effect  is  possible  because  the  radial  space  charge 
electric  field  as  defined  by  Gauss’  Law  can  be  easily  aligned  with  the  radial  electric  field  induced  by  the  generator’s 
magnet.  The  injected  space  charge  can  regulate  the  output  electric  field  while  the  magnetic  field  remains  constant. 
Space  charge  can  actually  reduce  the  amount  of  magnetic  field  needed,  thereby  reducing  the  magnet  current  or  even 
reducing  the  number  of  coil  windings  and  subsequence  magnet  volume  and  weight.  Space  charge  voltage  regulation 
does  not  require  significant  amounts  of  power  or  space  charge  current  due  to  the  low  fraction  of  the  overall  electron 
density.  The  space  charge  regulation  power  supply  is  only  a  fraction  of  the  mass  of  a  full  output  power  conditioner 
because  of  the  low  current  and  power  required  for  a  space  charge  power  supply.  This  reduces  the  energy  cost 
needed  to  fly  a  space  charge  power  supply  as  compared  to  the  energy  cost  needed  to  fly  a  full  output  power 
conditioner.  Large  technology  improvements  would  be  needed  to  increase  the  power  density  of  full  output  power 
conditioners  by  a  factor  of  5  or  10  to  make  them  light-weight  enough  to  not  serious  impact  the  weight  of  a  MHD 
power  system  (and  increase  its  mass  fraction  of  a  hypersonic  vehicle).  Replacing  traditional  external  power 
conditioners  with  an  internal  method  of  voltage  regulation  represents  a  major  technology  leap  in  MHD  generator 
technology  that  can  make  high  power  flight-weight  MHD  power  systems  practical  for  hypersonic  vehicles. 
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